Localization of non-interacting electrons in thin layered disordered systems 
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Localization of electronic states in disordered thin layered systems with b layers is studied within 
the Anderson model of localization using the transfer-matrix method and finite-size scaling of the 
inverse of the smallest Lyapunov exponent. The results support the one-parameter scaling hypothesis 
for disorder strengths W studied and fe = 1, . . . , 6. The obtained results for the localization length 
are in good agreement with both the analytical results of the self-consistent theory of localization and 
the numerical scaling studies of the two-dimensional Anderson model. The localization length near 
the band center grows exponentially with b for fixed W but no localization-delocalization transition 
takes place. 



I. INTRODUCTION 

The one-parameter scaling theory of Abrahams et ali 
predicts that all single-particle electronic states of non- 
interacting electronic systems in two dimensions (2D) are 
localized for arbitrarily weak random potential has at- 
tracted much attention in the last 25 years. This sub- 
tle effect is understood as being due to the construc- 
tive interference between time-reversed paths which in- 
creases the probability of returning to the original posi- 
tion leading to universal conductance fluctuations^ and 
corrections to conductivity, whose divergence in d = 2 
explains the absence of diffusion^iSiSiL2iS. It has also 
been investigated numerically for different lattice struc- 
tures, e.g., square, honeycomb and triangular at the band 
center {E — 0) and away from the band center in two- 
dimensional (2D) systems, but no extended states were 
found for the non-interacting and zero-field case, sup- 
porting the scaling hypothesises. There have been more 
theoretical and experimental investigations to see the va- 
lidity of the conclusions of scaling theorySsiiiiSiiii^iiS, 
supporting the theory from a variety of numerical, ex- 
perimental and analytical approaches in the absence of 
electron-electron interactions, magnetic field and spin- 
orbit interactiona^i'^1^1^1^1^^1^^. 

In 3D systems on the other hand, scaling theory pre- 
dicts that there exists a second-order continuous phase 
transition, with a critical disorder above which all the 
states are localized whereas below it extended states ap- 
pear at the band center, separated by mobility edges 
from the localized states in the rest of the spectrum. 
This was supported by numerical2i24iLi& and analyti- 
cal investigations'^^, with d = 2 being the lower criti- 
cal dimension of the transition. The question then can 
be posed what happens in thin films of finite thickness. 
When such a small thickness in terms of a few layers is 
introduced in an infinite 2D electron system, the degree 
of freedom of the electron motion will be increased. This 
thin film can be expected to behave as a 2D system, if 
the Thouless length Lxh, which is the length scale up to 
which an electron can diffuse without inelastic collisions, 
is larger than the film thickness 6, whereas for usual thin 
films, Lxh < 6 f^. The study of such systems should 



be important because experiments are always done in 
thin films of finite thickness, and some experimentsiS 
which studied electron transport with the variation of 
film thickness showed a fascinating role of the film thick- 
ness in the metal- insulator transition (MIT): the authors 
observed an MIT in Mo-C films at a certain thickness 
and found that the characteristic parameters like local- 
ization length and dielectric constant grow exponentially 
with increasing thickness near the transition. 

The numerical investigations in 2D and 3D sys- 
tems played an important role to study the proper- 
ties of electronic eigenstates in these systems. There 
have been various numerical methods applied in these 
systems, e.g. using exact diagonalization of the secu- 
lar matrices^, energy-level statisticsSi, transfer-matrix 
method (TMM)S4i2iii, studies of the dimensionless 
conduct ance22i24, etc. The results of these studies sup- 
port the one-parameter scaling hypothesis. 

Recently^^, there was a discrepancy between numer- 
ical results obtained using the TMM and an analytical 
analysis of thin films of finite thickness based on the self- 
consistent theory'*'5i. Numerically, an MIT induced by 
film thickness was found whereas the analytical results 
showed no metallic phase if the thickness of the film is 
finiteS^^. The discrepancy is rather surprising since previ- 
ous studies of 2D systems^i^iiSiS^ showed good agreement 
between the analytical and numerical results. 

With this motivation we address numerically the prob- 
lem of localization in thin disordered systems of nonin- 
teracting electrons using the standard TMM and finite 
size scaling (FSS)^'^^ for non-interacting electrons in zero 
magnetic field as described by the Anderson model of 
localization (defined below), in systems of considerably 
larger sizes than in Ref. |25l | and compare it with the 
analytical results. To this end we first note that the an- 
alytical expression for the localization length ^ from 
can be expressed as 



sinh b' 



jexp(6' C/W''), 



(1) 



where C, F are non-universal constants, W is the disorder 
strength and b' is the film thickness. The expression is 
valid for Ai? ^ ^ <C i < ixh and Xp b' L, where 



2 



Xp is the Fermi wavelength, £ the elastic mean- free path, 
and L is the system size. In infinitely thick films (6' — > 
cxd), on the other hand, a crossover from 2D to 3D takes 
place where one can get both a metallic and an insulating 
phase depending on the disorder strength. 

In section II, we describe the model and the notations; 
the TMM and FSS are discussed in section III; obtained 
numerical results are presented in section IV and some 
conclusions based on our numerical results drawn in sec- 
tion V. 



II. THE MODEL 

The Hamiltonian we have used for our numerical cal- 
culations is the Anderson model of localization given by 

N N 

H = Y,er\r){r\-tJ2 \r){r'\, (2) 

r {r.r') 

where |r) is a tight-binding state at site r. The sites r 
form a simple cubic lattice with N = M x b x L atoms 
where L is the length, M is the width, and b is the thick- 
ness of the lattice, here measured in the number of layers 
of the system, in units of the distance a between the 
layers. Since the only relevant microscopic scale in the 
continuum model that gives Eq. JQ) is i, we compare the 
length scales of the continuum and Anderson model via 
the relation b' — ab, where a = a/£ is a non- universal 
number which needs to be taken into account when com- 
paring analytical and numerical results. 

We use open boundary conditions along the b (thick- 
ness) direction and periodic boundary conditions along 
the M (width) direction. The on-site potentials 
are taken to be randomly distributed in the interval 
[—W/2,W/2]; {r,r') denotes that the hopping t is re- 
stricted to the nearest neighbors and the energy units 
are chosen by setting t = 1. 

We assume that this model describes the thin layered 
system if 6 ^ A/. When b becomes comparable to M at 
small W, one should expect a crossover from insulating 
towards metallic behavior after FSS, and it is therefore 
important to have b/M as small as possible in order to 
stay in the limit of finite thickness. We therefore at- 
tribute the appearance of the metallic-like behavior in 
the numerical simulations of Ref. 25] to the large val- 
ues of b/M used. Indeed, for b — 6, for instance, as M 
was changed from 2 to 15, the cross-section of the bars 
changed from rectangular to square to rectangular shape, 
and the scaling analysis for such geometries should not 
be compared with the analytical results of |2^. Instead, 
we focus on systems ^ with M ^ b. In the case of 
the thickest system studied here, for example, the scal- 
ing analysis is done for 6 = 6 and M — 23,33,47,65 as 
we discuss in detail below. 



III. ONE-PARAMETER FINITE SIZE SCALING 

In the TMM one determines localization properties by 
a scaling analysis of the Lyapunov exponents (LE) of the 
transfer matrices of quasi- ID bars. 

The largest length scale XMb is obtained from the 
smallest LE, XMb = ^/iMb- The accuracy of these Ajv/b 
is determined from the variance of the changes of the ex- 
ponents in the course of the iteration, and the relative 
error of A is 1% in our calculations. For 6=1 our system 
is a 2D strip of width M and length L ^ M, while for 
b — M , it becomes a 3D bar with square cross-section of 
M X M atoms and length L. 

We check the one-parameter scaling theory by the 
FSS analysis^ of the renormalized decay length, Am^ = 
XMb/M, for each b individually by expressing it as 

Ab{M,W)^hiUW)/M), (3) 

where ^b(W) is the scaling parameter which can be iden- 
tified as the localization length in the insulating regime 
and the resistivity in the metallic regime of the infinite 
system^ and fb{x) should be, if the one-parameter hy- 
pothesis holds, a universal function f{x) independent of 
b. The scale of ^ is determined from the condition that 
fb{x) = X when a; — > |8|, which corresponds to express- 
ing ^ in units of a. The numerical determination of fbix) 
is done by the method based on Ref. and the thus 
obtained scaling functions for different b are compared 
to check for the universality. If there is an MIT, then 
Af) for given b will have two branches, one which goes to 
zero as AI increases (indicating localization of electronic 
eigenstates) and the other which diverges as M increases 
(indicating extended eigenstates). 

Fig. n shows Aii{M,W) as a function of strip width 
M for various W, the scaling function fb{^b{W)/M), as 
well as the localization length ^{W) for systems with 
6 = 2, 3,4 layers, and Fig. |21 shows the same quantities 
for systems with & = 5, 6 layers. Our results for 6 = 1 
are not shown since they are in agreement with results 
obtained beforeSS. For larger W values, AMb decreases 
as the width of the system M increases indicating local- 
ization, while for smaller W values the decrease becomes 
much weaker. 

In addition to the results presented in Figs. ^ and 13 
we have calculated Af, for W values as small as 2, 2, 4, 4, 4 
for 6 = 2, 3, 4, 5, 6, respectively, but no metallic behavior 
was found. Instead, At{M) become approximately con- 
stant for smaller W values, which is in qualitative agree- 
ment with 2D being the lower critical dimension of the 
localization-delocalization transition. 

To see the dependence of the scaling function on b, 
we first plot all the obtained results for /f, in Figs. and 
|21using the same range on the axes for easier comparison, 
and then show all of the obtained ft together in Fig. |31 
without any rescaling or shifting of individual curves, in- 
cluding the curve for b = 1. The figure shows that fb 
does not depend on 6, therefore supporting the univer- 
sality of the scaling function f{x). Since it has only one 
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FIG. 1: Plot of the renormalizod decay length Kb(M,W) (left panels), scaling function fb (right panels) for systems of width 
M and 6 = 2,3,4 layers, from top to bottom. Dotted lines connect points for W values indicated in the corresponding legend. 
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FIG. 2: Same as Fig. Q for 6 = 5 and 6 layers, respectively. 
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FIG. 3: Plot of the scaling functions M^/M) for b = 
1,2,..., 6 from Figs, d and 



branch corresponding to the localized states it is also 
in agreement with the one-parameter scaling hypothe- 
sis for disorder strengths W > 4.75, 4.0, 4.5, 6.75, 7.25 for 
b — 2,3,4,5,6, respectively. For 6=1, one parameter- 
scaling was found®'^^ to hold for W > 2, and the pre- 
viously obtained^ scaling function f{x) for 2D system 
agrees within error bars with the f{x) obtained here for 
& > 1, furthermore supporting the universality of the re- 
sult. 

To compare the dependence of the localization lengths 
on W to the analytical result (0, we fit the obtained 
^b{W) to the form 



(4) 



as shown in the left panel of Fig. 01 The upper and lower 
right panels of the same figure show, respectively, the de- 
pendence of Cb and = Ab-^/sinh b/b on the number of 
layers b. The upper right panel of the figure also con- 



tains the mean least-squares linear fit of Cb, which gives 
Cb = C 6 Co, C = 99 ± 4, Co = 10 ± 15. The obtained 
value Co ~ is an additional check that the numerical 
results are in agreement with Eq. as far as the W~'^ 
dependence in the exponent is concerned. 

The obtained values of Ini^ depend on b, which seems 
to be in contradiction with Eq. but can be ex- 
plained by taking into the account the scale factor a, 
and Fig. |S1 shows ln_F, where F = Ai, ^ sinh ah jab for 
a = 0.73, when Ini^ becomes approximately constant, 
and a — 0.73 ± 0.13, where the error bar is estimated 
from the range of a values for which the values of In F for 
different b remain within error bars among themselves. 
This result can be interpreted as that the effective I of 
the lattice model ^ in the range of W studied is of the 
order of 1.4 ± 0.2 lattice constants. 

As an additional check of the validity of the above re- 
sults we notice that the obtained value of C is in agree- 
ment with the value C w 100 obtained in the numerical 
studies of the 2D Anderson mode&2S, i.e. for 6 = 1, even 
though the method by which we determine C relies es- 
sentially on the study of systems with 6 > 1. 



IV. CONCLUSIONS 

Thin quasi-lD bars with rectangular cross section of 
M X b and length L :s> M atoms were studied using the 
TMM and FSS of the inverse of the smallest Lyapunov 
exponent, for b — 1,2,3,4,5,6 and M ^ b. The ob- 
tained localization lengths were found to be in a good 
agreement with the analytical results. Despite the expo- 
nential increase of the localization length as the thick- 
ness increases, we find no signs of extended states for the 
disorder strengths studied, in agreement with the one- 
parameter scaling hypothesis. 
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